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A.  Introduction  and  Summary 

KMS  Fusion,  Inc.  began  investigating  the  two-stage  free  electron  laser 
(FEL)  concept  In  June  of  1980.  Since  then  we  have  developed  computer  simula¬ 
tion  codes  to  model  the  interactions  taking  place  between  the  electron  beam 
and  the  magnetic  and  electromagnetic  fields  in  the  two-stage  FEL  and  have 
Identified  parameter  regimes  In  which  such  a  device  Is  likely  to  operate 
successfully.  The  work  reported  here  is  In  support  of  a  two-stage  FEL  experi¬ 
ment  that  will  be  performed  at  the  University  of  California  Santa  Barbara 
(UCSB)  In  collaboration  with  Lius  Ellas  and  the  FEL  group  at  UCSB  using  their 
electrostatic  accelerator. 

In  a  two-stage  free  electron  laser  a  low  energy  (  ~  few  MeV)  electron 
beam  Is  used  to  produce  short  wavelength  (~  100-1000  M*n)  radiation.  This  Is 
done  In  two  steps.  First,  the  electron  beam  Is  passed  through  a  wlggler  mag¬ 
net  to  produce  long  wavelength  (~  100-1000  \m)  radiation.  This  radiation  Is 
then  backscattered  from  the  electron  beam  to  produce  the  short  wavelength 

radiation.  Useful  amounts  of  short  wavelength  radiation  can  only  be  produced 
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If  the  long  wavelength  electromagnetic  field  is  very  Intense  (~  10  W/cm  )  . 

During  the  past  year  we  have  developed  a  resonant  cavity  design  for 
containing  the  electromagnetic  pump  field  radiation,  which  provides  a  long, 
narrow,  high-intensity  region  for  the  FEL  Interaction  to  take  place  while 
expanding  the  beam  rapidly  outside  the  Interaction  region  to  minimize  the 
overall  cavity  length  needed  to  protect  the  end  mirrors.  This  quasi  optical 
cavity  design  Is  shown  In  Figure  1. 

The  design  has  a  number  of  unique  features.  A  low  loss  TE01  mode  is 
preferentially  propagated  In  the  cylindrical  waveguide.  This  mode  has  an 
annular  Intensity  distribution  with  intensity  minima  at  the  walls  and  on 
axis.  Holes  In  the  centers  of  the  cavity  end  mirrors  permit  the  second-stage 
short-wavelength  radiation  to  leave  the  cavity  with  negligible  pump  field 
losses  In  the  preferred  TEoi  mode* 

Because  the  TEqj,  mode  has  only  azimuthal  wall  currents.  It  can  propagate 
In  a  segmented  waveguide  structure.  Segmenting  the  waveguide  prevents  all 
modes  except  TEqM  modes  from  propagating  In  the  waveguide  and  also  permits 
Introduction  of  an  axial  electric  field  to  optimize  gain  In  both  the  first  and 
second  stages  of  the  FEL. 


Figure 


The  polarization  of  the  TE01  mode  Is  linear  but  varies  as  a  function  of 
axial  position.  A  helical  wlggler  magnet  Is,  therefore,  needed  to  fully 
excite  the  TE01  pump  field  mode.  For  best  coupling  to  the  pump  field  mode  It 
Is  also  desirable  to  use  an  annular  electron  beam. 

During  the  period  January  -  September,  1982  we  have  Investigated  this 
quasi optical  cavity  design  In  more  detail  and  have  begun  the  process  of  hard¬ 
ware  design.  We  have  also  started  developing  facilities  to  test  this  hardware 
and  characterize  It  before  It  Is  used  In  an  FEL  experiment.  In  addition,  we 
have  Investigated  problems  associated  with  electrostatic  accelerators  used  to 
drive  a  two-stage  FEL  and  have  Identified  accelerator  concepts  that  could  be 
scaled  to  high  average  currents  and  high  laser  power. 

In  this  report  we  summarize  the  work  that  has  been  completed  over  the 
past  nine-month  period.  In  Section  B,  laser  output  power  Is  calculated  for 
the  second  stage  of  a  two-stage  FEL  over  the  range  of  parameters  expected  In 
the  experiment  to  be  conducted  at  UCSB.  For  output  powers  up  to  1  kW,  second 
stage  output  power  was  found  to  Increase  as  the  cube  of  the  pump  field  Inten¬ 
sity  for  constant  electron  current,  and  as  the  sixth  power  of  the  electron 
current  when  the  pump  field  intensity  Is  assumed  to  increase  linearly  with 
electron  current. 

In  Section  C,  absorption,  diffraction,  and  mode  conversion  losses  for  the 
long  wavelength  pump  field  radiation  are  calculated  for  a  simple  quasi optical 
cavity  configuration.  Absorption  losses  were  found  to  be  about  0.51  per  round 
trip  and  mode  conversion  losses  about  1.6%.  Diffraction  losses  could  be  made 
less  than  0.2%.  The  cavity  considered  consisted  of  a  cylindrical  pipe  and  two 
mirrors.  We  believe  that  with  properly  designed  horns  at  the  ends  of  the 
waveguide  to  optimally  match  the  waveguide  to  free  space  modes,  mode  conver¬ 
sion  losses  could  be  greatly  reduced.  Reduction  of  absorption  losses  by 
cooling  the  cavity  structure  could  be  considered  If  mode  conversion  and  dif¬ 
fraction  losses  can  be  made  small  compared  with  absorption  losses.  To  measure 
cavity  loses  experimentally  before  Incorporating  the  cavity  Into  an  FEL  exper¬ 
iment,  we  are  setting  up  a  microwave  cold  test  facility.  This  facility  Is 
described  In  Section  D. 

For  high  efficiency,  continuous  FEL  operation,  the  wlggler  magnet  used  in 
the  first  stage  of  the  two-stage  FEL  would  have  to  be  either  superconducting 
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or  made  of  permanent  magnets.  Linear  higglers  are  commonly  made  using 
permanent  magnets,  but  helical  wlggler  magnets  employing  permanent  magnets 
have  not  yet  been  built.  We  are,  therefore.  Investigating  the  design  of  a 
permanent  magnet  helical  wlggler  structure.  Such  a  structure  could  be 
assembled  using  an  array  of  permanent  magnet  dipole  rings.  In  Section  E,  we 
describe  such  a  device,  determine  the  fields  that  could  be  produced  using 
SmCos  permanent  magnets,  and  Indicate  a  process  that  could  be  used  to  tune 
such  a  wlggler  to  correct  for  field  nonuniformities. 

In  order  to  perform  a  successful  two-stage  FEL  experiment,  several  defi¬ 
ciencies  of  the  present  UCSB  accelerator  must  be  corrected.  These  Include 
excessive  droop  and  ripple  In  the  accelerator  voltage,  an  Inadequate  power 
system,  and  the  ability  to  accept  only  a  very  small  (~  1%)  energy  spread.  In 
Section  F,  we  Indicate  steps  that  could  be  taken  to  correct  these  deficiencies 
so  that  a  proof  of  principle  two-stage  experiment  could  be  carried  out  at  the 
UCSB  facility.  We  also  present  a  design  concept  for  a  second  generation  elec¬ 
trostatic  accelerator  that  could  be  used  to  operate  a  high-power  two-stage  FEL. 

To  Improve  our  capability  to  model  FEL  systems  In  general  and  the  two- 
stage  FEL  experiment  In  particular,  we  have  developed  a  technique  for  modeling 
the  propagation  of  laser  radiation  in  a  free  electron  laser  In  two  and  three 
dimensions.  In  Section  G  of  this  report,  we  derive  a  three-dimensional  wave 
equation  for  the  electromagnetic  wave  in  an  FEL  using  a  basis  vector  formal¬ 
ism.  We  then  show  how  this  wave  equation  can  be  solved  using  Fourier  trans¬ 
forms  In  Cartesian  coordinates.  The  technique  is  extended  to  cylindrical 
coordinates,  but  problems  arise  because  fast  Fourier  transform  techniques  can¬ 
not  be  used  to  obtain  an  Incremental  solution  to  the  radial  part  of  the  wave 
equation.  An  approximation  that  permits  the  use  of  the  fast  Fourier  trans¬ 
forms  Is  Introduced  and  a  method  Is  developed  to  Insure  the  accuracy  of  the 
solution  over  a  finite  spatial  mesh.  A  simplified  two-dimensional  model  of 
the  Interaction  of  the  electron  beam  with  an  axl symmetric  radiation  field 
profile  Is  used  to  obtain  Information  on  the  effect  of  radial  variation  of  the 
laser  field  on  trapping. 

B.  Two-Stage  FEL  Scaling 

Basic  scaling  relationships  for  a  two-stage  FEL  have  been  obtained  using 
a  simple  one-dimensional  resonant  particle  amplifier  simulation  code.  The 
experimental  variables  that  were  found  to  have  the  greatest  effect  on  second 


stage  output  power  are  the  intensity  of  the  electromagnetlc(em)  pump  field 
produced  In  the  first  stage,  the  electron  beam  current,  and  second-stage 
optical  system  losses.  We  have  calculated  laser  gain  over  a  range  of  values 
for  the  pump  field  Intensity,  laser  field  Intensity,  and  electron  current 
assuming  the  resonant  cavity  structure  of  Figure  1  and  using  accelerator  and 
resonant  cavity  parameters  that  are  expected  for  the  UCSB  two-stage  FEL 
experiment. 

In  our  calculations  we  assume  the  electron  energy  is  3  MeV,  the  frac¬ 
tional  electron  energy  spread  Is  10"4  ,  and  elec  jn  current  Is  In  the  range 
between  2  and  20  amp.  The  length  of  the  second  stage  interaction  region  Is 
assumed  to  be  3  m  and  the  pump  field  wavelength  Is  assumed  to  be  1  mm,  which 
gives  a  laser  wavelength  of  5.3  um.  The  diameter  of  the  cylindrical  waveguide 
Is  assumed  to  be  2.44  cm. 

An  annular  electron  beam  Is  assumed  to  Interact  with  the  region  of  peak 
Intensity  of  the  em  pump  field  (see  Figure  2).  The  cross  sectional  area  of 
the  Interaction  region  Is  assumed  to  be  1.26  cm2.  In  Figure  3,  second  stage 
output  power,  Pout,  Is  plotted  as  a  function  of  laser  power  at  the  input  end 
of  the  second  stage  Interaction  region,  Pun,  f°r  a  number  of  values  of  the 
pump  field  power.  Pout  Is  obtained  assuming  all  of  the  laser  light  produced 
In  the  amplifier  Is  available  as  output  power.  This  Is  probably  a  reasonably 
good  assumption  If  several  percent  of  the  light  Is  removed  from  the  cavity  on 
every  pass.  If  total  cavity  losses  must  be  kept  below  a  percent,  then  a 
significant  fraction  of  the  light  produced  may  not  be  available  as  useful 
output  power. 

In  Figure  3  we  see  that  laser  output  power  Increases  both  with  Increasing 
pump  field  Intensity  and  increasing  input  laser  intensity.  In  an  FEL  oscilla¬ 
tor  at  Its  equilibrium  operating  point,  at  which  optical  system  losses  just 
equal  laser  gain,  the  Input  laser  power  and  pump  field  power  will  not  be  inde¬ 
pendent  quantities.  Superimposed  on  Figure  3  are  two  curves  for  1%  and  10% 
amplification.  If,  the  optical  system  is  designed  to  produce  10%  loss  per 
pass,  only  the  region  of  parameter  space  above  the  10%  gain  curve  will  be 
accessible.  This  places  a  lower  limit  on  the  pump  field  power  needed  to 
obtain  a  given  laser  output  power.  For  example,  with  the  parameters  used  In 
Figure  3,  If  the  pump  field  power  Is  108  watts  and  optical  system  losses  are 
10%  per  pass,  the  maximum  output  power  that  could  be  attained  would  be  0.45  watt. 


E(r)  is  within  10%  of  Emax  in  interval  .34  <  r/a  <  .62 
Bucket  height  (proportional  to  E(r)fy  varies  by  5%  over  same  interval 
Cross  sectional  area  of  electron  beam  •  .27  of  waveguide  cross  section 
<E{r)>  ■  .96  Emax  in  interaction  region.  E^  occurs  at  r/a  -  .47 
Example:  a  ■  1 .22  cm,  guide  area  *  4.68  cm3 
electron  beam  area  ■  1.26  cm3 


Figure  2 


OUTPUT  POWER  SCALING  FROM  EM  PUMP  FEL  (UCSB  PARAMETERS) 


Figure 


In  Figure  4  output  power  Is  plotted  as  a  function  of  pump  field  power  for 
the  parameters  of  Figure  3.  Two  curves  are  shown,  one  for  10%  loss  per  pass 
in  the  second-stage  optical  system,  the  other  for  1%  loss.  The  slopes  of 
these  curves  Indicate  that  laser  output  power  scales  as  the  cube  of  the  pump 

field  power.  We  also  see  that  with  a  two-amp  beam  very  Intense  pump  fields 

are  needed  to  obtain  even  modest  laser  output  powers. 

Pump  field  power  is  determined  by  dividing  the  power  gain  of  the  pump 

field  by  the  fraction  of  the  energy  lost  per  cavity  round  trip.  Input  power 

to  the  pump  field  Is  limited  by  generator  capacity  In  the  accelerator  and  the 
maximum  energy  with  which  an  electron  can  lose  and  still  return  to  the  dome  of 
the  accelerator. 

The  generator  In  the  accelerator  at  the  present  time  provides  20  kw  of 
power  In  the  dome.  This  Is  sufficient  to  provide  for  an  average  of  10  keV 
energy  loss  per  electron  for  a  2-amp  beam  If  all  of  the  power  can  be  used  for 
this  purpose.  Electron  gun  voltage  Is  50  kV  so  that  losses  of  up  to  a  few 
tens  of  keV  for  some  electrons,  with  an  average  10  keV  /  electron  loss  for  the 
entire  beam.  Is  acceptable  for  collection  of  the  return  beam. 

Round  trip  cavity  losses  are  Initially  expected  to  be  between  1  and  .1% 
for  the  pump  field  radiation.  Losses  may  be  reduced  below  .1%  at  a  later  time 
with  Improvements  In  cavity  design.  With  20  kw  of  Input  power  and  cavity 
losses  of  .1%  per  pass,  the  Intracavity  power  would  be  2  x  1C7  watts.  In 
order  to  obtain  pump  field  powers  of  108  watts  or  greater,  more  generator 
power  must  be  provided  for  the  electron  beam,  a  greater  energy  loss  per  elec¬ 
tron  must  be  permitted,  and  cavity  losses  must  be  reduced.  It  Is  expected 
that  all  of  these  things  will  happen  as  operating  experience  Is  gained  an  the 
system  Is  upgraded. 

It  Is  presently  planned  to  modify  the  UCS8  accelerator  to  Increase  the 
maximum  beam  current  from  2  to  20  amp.  It  Is,  therefore.  Important  to  know 
how  laser  output  Is  expected  to  scale  with  increasing  beam  current.  Several 
parameters  will  change  as  a  result  of  Increasing  the  beam  current. 

Increasing  the  electron  beam  current  will  Increase  the  gain  In  the 
first  stage  of  the  laser  as  well  as  In  the  second  stage.  Therefore,  the 
pump  field  Intensity  will  be  greater  than  could  be  achieved  with  a  lower 
beam  current.  If  the  fraction  of  the  electron  energy  converted  to  long 
wavelength  radiation  remains  the  same,  the  pump  field  power  will  Increase 
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linearly  with  electron  current.  It  may  be  possible  to  increase  the  pump 
field  power  more  rapidly  than  linearly,  however,  since  conversion  effi¬ 
ciency  could  increase  at  the  higher  pump  field  intensity. 

For  the  pump  field  power  to  increase,  the  power  supplied  by  the 
generator  must  also  increase.  This  might  be  done  by  increasing  the 
generator  capacity,  storing  generated  power  between  pulses,  or  a  combi¬ 
nation  of  the  two.  In  determining  the  scaling  of  second-stage  output 
power  with  beam  current,  we  will  make  the  conservative  assumption  that 
the  pump  field  power  Increases  only  linearly  with  electron  current  due 
to  accelerator  system  limitations.  Accelerator  limitations  will  be  less 
of  a  factor  for  second  stage  gain  in  the  UCSB  experiment,  provided 
conversion  efficiency  In  the  second  stage  is  small  compared  to  the  first 
stage. 

The  second-stage  output  power  will  increase  much  more  rapidly  than 
linearly  with  Increasing  beam  current,  because  both  the  pump  field 
intensity  and  the  second  stage  gain  at  a  given  pump  field  intensity  will 
increase.  If  the  optical  system  has  a  fixed  percentage  loss  per  pass, 
the  equilibrium  operating  point  will  shift  upward  as  the  curve  of  Pout 
as  a  function  of  PL1n  rises  relative  to  a  line  of  constant  gain  (see 
Figure  3). 

The  effects  of  Increased  pump  intensity,  increased  second  stage 
gain  at  a  given  pump  field  Intensity,  and  shift  of  the  equilibrium 
operating  point  are  all  Included  in  the  curves  of  laser  output  power  as 
a  function  of  electron  current  plotted  In  Figure  5.  Curves  are  shown 
for  optical  systems  with  both  10%  and  1%  losses  per  pass.  For  peak 
output  powers  below  several  kilowatts,  Pout  scales  approximately  as  the 
sixth  power  of  the  electron  current.  Above  10  kW  output  power  Increases 
more  slowly  with  Increasing  current.  The  reason  for  the  difference  in 
scaling  Is  that  at  the  lower  laser  powers  the  fraction  of  the  electron 
distribution  contributing  to  laser  gain  increases  as  the  intracavity 
laser  Intensity  Increases.  Above  a  certain  power  level  the  phase  space 
buckets  are  large  enough  to  trap  most  of  the  electrons  In  the  beam,  so 
that  the  fraction  trapped  remains  constant. 

The  output  power  shown  In  Figure  5  Is  the  peak  power  that  could  be 
obtained  In  the  UCSB  experiment  for  the  assumed  resonant  cavity  losses. 
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SCALING  OF  OUTPUT  POWER  WITH  ELECTRON 
CURRENT  FOR  UCSB  TWO-STAGE  FEL 


Figure  5 


Since  the  accelerator  will  lose  charge  at  a  more  rapid  rate  than  it  can 
be  replenished  by  the  input  current,  the  accelerator  will  operate  in  a 
pulsed  mode.  During  each  electron  pulse,  both  the  pump  field  and  the 
laser  field  must  build  up  from  spontaneous  emission. 

If  the  time  needed  to  reach  equilibrium  is  longer  than  the  electron 
pulse,  the  peak  laser  output  power  will  be  less  than  the  equilibrium 
value.  Electron  pulse  length  is  directly  dependent  on  the  efficiency 
with  which  the  electron  beam  is  returned  to  the  dome  of  the  accelerator. 
If  IX  of  the  electron  beam  is  lost  in  the  beam  transport  system  the 
electron  pulse  length  would  be  about  80  psec  for  a  2-amp  beam  in  the 
UCSB  accelerator  and  about  8  psec  for  a  20-amp  beam.  Initial  experiments 
have  produced  collection  efficiencies  of  about  96%  without  an  PEL 
amplifier  in  the  beam  line.  Collection  of  99%  or  better  with  the  PEL  in 
operation  are  a  long  range  goal.  We  have  assumed  an  annular  electron 
beam  in  the  calculation.  The  first  experiments  will  not  be  done  with  an 
annular  beam,  and  therefore  only  part  of  the  electron  beam  will  interact 
with  the  high  intensity  pump  field,  lowering  the  effective  electron 
current. 

In  conclusion,  we  have  found  that  laser  output  power  should  vary 
as  the  sixth  power  of  the  electron  current  in  the  UCSB  two  stage  FEL 
experiment.  If  the  electron  current  is  held  constant,  laser  output 
power  will  vary  as  the  third  power  of  the  pump  field  Intensity. 

Decreasing  optical  system  losses  increases  the  amount  of  laser  light 
produced  in  the  cavity.  Optical  system  losses  include  both  the  output 
power  and  dissipative  losses  such  as  mirror  absorption  and  scattering. 
Dissipative  losses  In  the  second  stage  can  probably  be  kept  to  a  percent 
or  less  per  round  trip  optical  pass. 

The  scaling  with  pump  field  power  and  electron  current  for  the 
two-stage  FEL  Indicates  that  prospects  for  developing  a  high-power 
hlgh-efflclency  device  based  on  this  concept  appear  very  good.  If  the 
UCSB  experiment  works  as  predicted,  high  average  power  devices  could  be 
built  with  electron  currents  not  exceeding  20  amp.  This  would  require 
developing  an  accelerator  system  with  sufficiently  low  losses  and  suffi¬ 
ciently  high  Input  current  to  operate  In  a  D.C.  mode.  It  would  also 
require  much  higher  power  conversion  from  the  electron  beam  to  the 


optical  beam  and  higher  pump  field  power  levels.  As  will  be  discussed 
in  the  following  sections,  this  scaleup  may  be  accomplished  with  existing 
state  of  the  art  technology. 


C.  Resonant  Cavity  Design 

The  pump  field  cavity  must  be  designed  to  minimize  cavity  losses. 
Cavity  losses  reduce  the  maximum  pump  field  intensity  that  can  be 
attained  in  the  interaction  region  and  also  reduce  the  overall  effi¬ 
ciency  of  the  device.  There  are  three  main  classes  of  losses  that  can 
occur  in  the  cavity  —  absorption  losses,  diffraction  losses  and  mode 
conversion  losses. 

1.  Absorption  Losses:  For  good  conductors,  the  surface  resistance 
Rc  can  be  written  as1 


Rs  58  06 


(C-l) 


where  o  is  the  conductivity  and  6  is  the  skin  depth.  The  skin  depth  is 
gi ven  by 


6  *  (*cn0f) 


(C-2) 


where  f  Is  the  frequency  of  the  wave.  For  normal  Incidence  illumination 
It  can  be  shown  theoretically2’ 3  and  experimentally k  that  the  absorption 
loss  Is  given  by 

Pr  »  P1  (1  -  4  R$/Z0)  (C-3) 


where  Zo  »  Po/£o  *  3770  is  the  free  space  Impedance,  Pr  Is  the  reflected 
power,  and  Pj  is  the  Incident  power. 

Table  1  gives  conductivity,  skin  depth,  surface  resistance,  and  the 
dimensionless  quantity  R$/Z0  for  a  number  of  different  metals  at  a 
wavelength  of  1  mm.  The  values  In  Table  I  are  based  on  the  assumption 
of  idealized  surface  conditions.  Surface  contamination  and  Imperfections 
on  grain  boundaries  could  Increase  absorption  by  up  to  a  factor  of  2. 

For  the  case  of  a  cavity  with  two  copper  mirrors  at  room  temperature 
the  round  trip  cavity  transmission  factor  Tm  is  given  by 


Tm  •  1  -  8  Rs/Z0  -  0.997 


(C-4) 
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Table  I 


NORMAL  INCIDENCE  ABSORPTION  LOSSES  ON  THE  SURFACE 
OF  GOOD  CONDUCTORS  AT  ROOM  TEMPERATURE 

(X  *  1  mm) 


Conductor 

Conductivity 
a  (mho/m) 

Skin  depth 

6  (m) 

Surface 

Resistance 

Rs  (0)  _ 

Loss 

Factor 

Rs/zo 

A1 uml num 

3.54  x  107 

1.54  x  10"7 

0.18 

4.9  x  10‘4 

Copper 

5.8  x  107 

1.21  x  10*7 

0.14 

3.79  x  10"** 

SI 1 ver 

6.15  x  107 

1.17  x  10-7 

0.139 

3.68  x  10"4 

Gold 

4.5  x  108 

1.37  x  lO"7 

0.16 

4.3  x  10~4 
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In  a  cylindrical  waveguide  the  attenuation  of  the  power  in  a  propa¬ 
gating  wave  is  given  by 


P(z) 


Pi  e 


-ca 


(C-5) 


For  a  TE01  wave  the  attenuation  coefficient,  a,  is  given  by2*5 


2R 


(fc/f): 


aZ0  v'l  -  ( fc/f )2 

where  a  is  the  waveguide  radius,  f  is  the  frequency  of  the  radiation, 
and  fc  is  the  waveguide  cutoff  frequency.  For  a  TEqi  mode 


(C-6 ) 


u  X. 
2n  a 


(C-7) 


where  u  is  the  value  of  the  first  off-axis  zero  of  the  Jx  Bessel  function, 
Ji(u)  *  0  for  u  *  3.83.  For  X  =  .1  cm  and  a  *  1.22  cm,  fc/f  *  0.05 
and  a  *  1.55  *  10“**  m“l.  For  a  5-m-long  copper  waveguide,  at  room 
temperature,  the  round  trip  transmission  factor  Tg  will  be  given  by 


-2aL 

e 

.998 


(C-8) 


Total  transmission  losses  due  to  absorption  In  both  the  mirrors  and  the 
waveguide  in  one  cavity  round  trip  will  therefore  be 


Absorption  loss  ■  1  -  TmTg  *  1-.995  *  5xl0"3  (C-9) 

Whether  a  loss  rate  this  high  would  be  tolerable  in  a  high  power 
two-stage  FEL  depends  on  the  fraction  of  the  electron  energy  converted 
to  usable  photon  energy  per  pass,  the  overall  laser  efficiency  required, 
and  the  size  of  the  losses  compared  with  other  system  losses.  Absorption 
losses  could  be  made  almost  arbitrarily  low,  if  required,  by  cooling  the 
metallic  surfaces  of  the  cavity.  Cooling  should  only  be  considered 
after  It  is  demonstrated  that  absorption  losses  are  the  dominant  cavity 
loss  mechanism. 


2.  Mirror  Diffraction  Losses:  If  the  radius  of  a  free  space  mode 
leaving  the  waveguide  is  w0,  then  the  radius  of  that  mode  at  any  distance 
z  from  the  waveguide  is  given  by 


U)(z)  =  wo  (  1  +  (yjjO  )  ,Z 


(C-10) 


where 


Z0  »  — r (C-ll) 

^  ★ 

If  we  assume  the  intensity  pattern  is  primarily  that  of  a  TEMoi  annular 
mode,  the  intensity  across  the  beam  is  given  by 


I ( r )  *  1(0)  pe 


(C-I2) 


where 


(C-13) 


The  value  chosen  for  wo  Is  some  fraction  *  of  the  waveguide  radius,  a, 
which  gives  a  good  match  between  the  free  space  and  waveguide  cavity 


modes. 


(C-H) 


If  we  assume  that  the  TE0i  waveguide  mode  Is  converted  into  a  TEMoi 
free  space  mode  then  a  best  fit  value  for  o  Is  a  «  0.5.  For  a  »  1.22  cm, 
X  »  0.1  cm,  and  an  assumed  distance  of  2  m  between  the  end  of  the 
waveguide  and  the  mirror,  the  beam  radius  at  the  mirror,  w^,  will  be 
10.45  cm. 


/  pe’p  dp  -  -(1+p)  e" 


(C-15) 


and  the  value  of  this  Integral  from  0  to  *  equals  1.  Therefore  to  find 
the  fraction  of  the  cavity  radlaton  Inside  a  mirror  of  radius  A  we  only 
have  to  evaluate  (C-15)  from  0  to  p(A).  For  example,  If  we  chose  A  »  2u^, 
then  p(A)  »  8  and  .997  of  the  radiation  will  hit  the  mirror.  If  A 
»  2.5  then  .99995  of  the  radiation  leaving  the  waveguide  will  Intercept 
the  mirror.  Therefore,  diffraction  losses  around  the  outside 


of  the  mirror  can  be  made  arbitrarily  small  by  making  the  mirror  suffi¬ 
ciently  large. 

The  short  wavelength  radiation  will  also  come  out  of  the  waveguide  as 
an  annular  beam,  but  it  will  spread  much  less.  For  5.3  pm  radiation  the 
beam  radius  at  a  mirror  2  m  from  the  waveguide  would  be  0.61  cm, 
essentially  the  same  as  at  the  waveguide.  If  we  make  a  hole  in  the  mirror 
with  a  radius  2.5  times  larger,  virtually  no  short  wavelength  radiation 
will  hit  the  mirror.  From  (C-15)  the  fraction  of  the  long  wavelength 
radiation  that  would  be  lost  through  this  hole  would  be 
1  *  10"3  with  a  round  trip  loss  due  to  diffraction  of  2  x  10"3.  This  loss 
could  be  reduced  by  placing  the  mirrors  farther  from  the  waveguide.' 

3.  Mode  Conversion  Losses:  During  the  period  covered  by  this  report 
a  model  has  been  developed  to  calculate  losses  from  the  simplest  possible 
quasloptlcal  cavity  consisting  of  a  pipe  and  two  mirrors.  In  this  model 
a  TEoi  waveguide  mode  is  converted  to  a  linear  combination  of  TEMpi  modes 
at  the  end  of  the  waveguide.  The  individual  free  space  modes  propagate  to 
the  mirror  and  back.  The  returning  power  feeds  the  TEoi  mode  as  well  as 
other  waveguide  modes  and  is  partially  lost  around  the  outside  of  the 
waveguide.  Losses  are  assumed  to  consist  of  all  energy  not  reconverted 
into  the  TEoi  mode. 

The  discontinuity  at  the  end  of  the  waveguide  results  in  reflection 
of  part  of  the  field  back  Into  the  waveguide.  For  a  simple  pipe  the 
reflection  coefficient  Rg  is  given  by 

Rg  »  (Z0  -  Zg)  /  (Z0  +  Zg)  (C-16) 

where 

Zg/Z0  ■  (1  -  (fc/f)J)1/2  (C— 17) 

For  the  case  we  are  considering  Rg  *  6.25  x  10"4.  This  reflection  need 
not  degrade  cavity  performance  If  the  reflected  wave  is  in  phase  with  the 
wave  reflected  from  the  cavity  end  mirror.  We  will  assume  this  to  be  the 
case.  Improved  cavity  designs  with  conical  horns  at  the  end  of  the 
waveguide  could  significantly  reduce  this  reflection. 

Inside  the  waveguide  the  electric  field  has  the  form 


(C-18) 


£Q(r,z,t)  *  Eq{ r )  cos  (kz  -  wt) 
where  EQ(r)  Is  proportional  to  the  J1  Bessel  function,  ^(3.832  r/a). 


At  the  end  of  the  waveguide  the  Bessel  function  mode  Is  reexpressed 
as  a  linear  combination  of  Gauss  LaGuerre  modes  with  the  azimuthal 
component  of  the  electric  field  given  by 


E0p(u)  pl)(u> 

_L  e-u/2 
yp+r 

(C-19) 

where 

u  --§4 

a2  a2 

( C — 20 ) 

and 

a  *  w0/a 

(C-21) 

u)0  Is  the  beam-waist  radius  of  the  free  space  beam  at  the  end  of  the 
guide.  The  associated  LaGuerre  polynomials  are  obtained  from  the  relation 

L£l)(u)  *  ‘  if  Lp+l(u)  (C~22) 

The  waveguide  and  free  space  modes  are  normalized  so  that 

J*  2*r  dr  E2(r)  *  1  (C-23) 

0  B 
and 

/'  2«r  dr  E|p  ■  1  (C-24) 

At  the  boundary  the  waveguide  mode  can  be  expressed  In  terms  of  the 
free  space  modes  In  the  following  way: 

E9  *  p^O  CP  E®P  (C-25) 

where 

Cp  ■  J*  2 nr  dr  Ee(r)  Eep(2r2/a2a2)  (C-26) 

A  computer  program  has  been  written  to  solve  these  equations.  Cases  were 
considered  with  up  to  12  free  space  modes.  A  range  of  values  for  the 
parameter  *  were  tried.  For  12  modes  a  best  fit  was  obtained  for 


> 


18 


a  »  0.3.  The  values  of  the  coefficients  Cn  for  this  case  are  given 
in  Table  2.  The  sum  of  the  squares  of  the  first  12  terms,  Cp2,  is 
also  given.  The  deviation  of  this  sum  from  unity  is  a  measure  of  the 
contribution  due  to  modes  with  p  >  11.  Cr  is  the  remaining  amplitude 
if  only  one  higher  order  mode  were  needed  to  completely  represent  the 
waveguide  mode  at  the  boundary. 


Table  II 


Coefficients  for  Expansion  of  TEoi  modes  as  a  series  of 
Gauss  Laguerre  free  space  modes  for  a  *  0.3 


Co  * 

.6153 

C6  a  -.3733 

Cl  a 

-.5828 

C7  *  .2040 

C2  - 

.4400 

C8  -  .1125 

C3  ■ 

-.2695 

C9  a  -.1388 

Cl,  a 

.1139 

Cj^q  *  -.5760 

C5  * 

-.4748 

x  10- 2 

Cu  •  .8875 

11 

Cr  - 

.1473 

x  10-1 

l  C  2  *  .9998 

p 

p-0  p 

x  10'1 
x  10-1 
x  10”1 
x  10”1 
x  10"2 
x  10"  2 


When  the  radiation  leaving  the  guide  propagates  to  the  mirrors  and 
back  to  the  waveguide,  the  different  modes  do  not  stay  In  phase.  The 
phase  shift  of  the  modes  with  p>0  relative  to  the  p*0  mode  Is  given  by 

♦  »  4p  arctan  (  ^  -j)  (C-27) 

where  Is  the  distance  from  the  end  of  the  waveguide  to  the  mirror. 

For  Lgm  *  200  cm,  X  >  .1  cm,  a  ■  1.22  cm  and  a  *  .3 

♦_ 

■  arctan  (47.5)  ■  88.8  degrees  (C-28) 
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The  modes  with  small  values  of  p  win  have  relative  phase  shifts  close  to 
zero  (multiples  of  2n  radians)  while  modes  with  very  large  p  will  have 
small  amplitudes. 

The  radiation  returning  to  the  waveguide  is  given  by 


l  CnEflne 
p=0  p  9p 


j=i 


lij 

V  J»j 


n=0 


qnE0n 


(C-29) 


In  (C-29)  the  Incoming  wave  is  decomposed  Into  two  parts.  The  first  sum 
on  the  right  hand  side  are  the  waveguide  modes  which  may  be  nonzero  for 
r<a  but  are  equal  to  zero  for  r  >  a.  The  second  sum  on  the  RHS  of  (C-29) 
represents  the  remainder  of  the  radiation  field  outside  the  waveguide  and 
equals  zero  for  r  <  a. 

Multiplying  (C-29)  by  2ttrp|g  CpEep  and  Integrating  using  the  ortho¬ 
gonality  conditions  gives 


*  1$  1«i 

l  C2  e  P  -  ALe  (C-30) 

p»0  P 

This  result  Is  obtained  because  p|Q  CpEep  Is  equivalent  to  the  TE01  mode 
In  the  waveguide  (Bi-Eg),  and  all  the  other  waveguide  modes,  B-,-  (1>1),  are 
orthogonal  to  It.  Outside  the  waveguide  p|Q  CpEep*  0,  so  the  Integrand 
for  the  second  sum  In  (C-29)  Is  always  zero.  The  fraction  of  the 
returning  power  converted  back  Into  the  TEoi  mode  Is  equal  to  AL2. 


Ai2  -  [Re[Aiet®i]]2  +  [yAie^M] 


(C-32) 


(  I  C  2  cos*)2  +  (  l  C  2  sin*  ) 


•0  p 


0  P 


* 


Taking  only  the  first  12  terms  of  each  sum  we  get  a  value  for  AL2  «  0.9921. 
This  Indicates  that  practically  all  of  the  radiation  returning  to  the 
waveguide  will  reenter  the  guide  and  again  propagate  as  a  TE01  mode. 

Upper  and  lower  bounds  on  the  value  that  would  be  obtained  If  an  Infinite 
sum  were  taken  are  ■  AL2  ±  2Cr2.  The  round  trip  transmission  factor  for 
the  TE0i  mode  equals  AL4  since  mode  conversion  would  occur  at  both  ends 
of  the  waveguide.  For  the  first  12  modes  A]4  *  .9843  ±  .0008. 
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A  loss  of  .8%  per  round  trip  pass  due  to  mode  conversion  is  obvi¬ 
ously  higher  than  we  would  like,  but  this  result  is  only  the  first  cut 
at  modeling  an  initial  unoptimized  configuration.  Lower  losses  could 
be  obtained  with  this  same  configuration,  however,  by  increasing  the 
distance  between  the  end  of  the  waveguide  and  the  mirrors  with  a  corre¬ 
sponding  increase  in  the  size  of  the  mirrors. 

D.  Measurement  of  Cavity  Q 

Modeling  calculations  of  the  quasioptical  cavity  indicate  that 
losses  per  round  trip  pass  could  be  kept  below  1%.  We  wish  to  verify 
these  calculations  experimentally  before  fabricating  the  hardware  to  be 
used  in  the  UCSB  two-stage  experiment.  To  do  this  we  are  setting  up  a 
microwave  cold  test  laboratory  in  which  the  Q  of  cavities  similar  to 
that  which  will  be  used  in  the  two-stage  FEL  experiment  can  be  measured. 
The  testing  procedures  developed  in  this  laboratory  will  not  only  provide 
us  with  a  better  understanding  of  cavity  performance,  but  will  also 
provide  us  with  diagnostic  experience  that  will  be  needed  for  the  FEL 
experiment. 

We  have  determined  the  cold  test  measurements  which  will  be  neces¬ 
sary  during  the  construction  and  subsequent  characterization  of  the 
quasi -optical  pump  cavity.  These  consist  mainly  of  progressively 
narrower  bandwidth  measurements.  The  test  cavity  configuration  shown 
in  Figure  6  will  be  used  for  testing  cavity  components.  The  information 
needed  for  the  design  of  the  pump  cavity  can  be  obtained  from  the  test 
cavity. 

The  cold  testing  will  be  performed  at  a  wavelength  of  3  mm.  An 
efficient,  reliable  mode  converter  to  launch  the  low  loss  TE01  mode  into 
the  circular  guide  is  available  at  3  mm,  but  not  at  shorter  wavelengths. 
Also,  good  frequency  stability  (a  few  parts  in  108}  is  required  of  the 
bench  source.  This  stability  can  be  guaranteed  at  3  mm,  but  again  not 
at  shorter  wavelengths. 

Two  bench  sources  will  be  used  to  cover  the  range  of  bandwidths  we 
need.  A  frequency-swept  backward  wave  oscillator  will  provide  the  wide 
bandwidth  necessary  for  many  of  our  initial  tests  and  will  permit  the 
measurements  to  be  made  quickly.  These  tests  include  identification  of 
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CONVERTOR 


Figure 


undesirad  inodes  and  initial  Q  measurements  of  the  desired  mode.  Q  can 
be  expressed  as  f/Af,  where  f  is  the  center  frequency  of  the  mode,  and 
Af  is  its  bandwidth,  which  is  proportional  to  the  losses  suffered  by  the 
mode. 

The  test  cavity  is  expected  to  have  absorption  losses  of  a  few 
tenths  of  a  percent  per  round  trip  pass  at  room  temperature.  A  compa¬ 
rable  amount  of  power  will  be  transmitted  through  the  element  coupling 
the  cavity  to  the  cold  test  equipment.  This  coupler  will  be  designed 
so  as  not  to  produce  undesired  modes.  We  expect  that  losses  at  the 
waveguide/free-space  transition  will  be  less  than  0.1%  for  optimized 
cavities.  The  bandwidth  corresponding  to  the  total  losses  will  be 
several  parts  In  107.  Thus,  In  order  to  detect  improvements  In  the 
waveguide/free-space  transition,  we  will  need  frequency  stability  of 
a  few  parts  in  10s.  A  phase  locked  Gunn  diode  oscillator  meets  this 
requirement  and  has  an  adequate  power  level. 

An  alternative  method  of  measuring  cavity  Q  Is  by  measuring  the 
1/e  decay  time  (tc  *  Q/2itf)  of  the  natural  oscillations  of  the  cavity. 
This  may  be  done  with  a  high  isolation  switch  (shown  in  Figure  6)  and 
sensitive  diode  detector. 

E.  Permanent  Magnet  Helical  Wiggler  Design 

A  helical  wiggler  magnet  is  needed  to  fully  excite  the  TE0i  pump 
field  mode.  During  this  contract  period  we  have  investigated  the 
possibility  of  producing  this  helical  wiggler  field  using  an  array  of 
permanent  magnets.  We  have  found  that  it  is  possible  to  do  this  and 
that  the  magnetic  field  produced  using  permanent  magnets  would  have  the 
same  form  as  that  which  Is  produced  using  a  bi filar  helical  winding  such 
as  In  the  Stanford  wiggler6. 

1.  Multi  pole  Rings:  The  building  block  for  a  permanent  magnet 
(SmCos)  helical  wiggler  will  be  a  dipole  ring.  The  design  concept  for 
any  multipole  ring  is  as  follows.  If,  In  a  ring  of  perfect  rare  earth 
xobalt  (REC )  material  with  Inner  radius  rj,  and  outer  radius  r2,  the  easy 
axis  Is  continuously  rotated  according  to  the  relation 

a  -  (N+l)e  ,  (E-l) 


aan./.: 


then  a  perfect  multi  pole  field  of  harmonic  order  N  will  be  produced 
inside  r^7  The  strength  of  these  multipole  fields  is,  for  N  >  2, 


BN(ri)| 


(E-2) 


where  the  harmonic  number,  N,  equals  half  the  number  of  poles  in  the 
magnet  ring.  For  N  =  1,  the  dipole  field  strength  is 


(E-3) 


It  has  been  shown7  that  this  multipole  ring  design  produces  the  strongest 
possible  field  for  a  given  amount  of  REC  material.  It  is  not  presently 
possible  or  practical  to  produce  rings  with  a  continuous  easy  axis 
orientation.  The  solution  is  to  segment  the  ring  into  a  number  of 
pieces  as  illustrated  in  Figure  7.  In  each  segment  the  easy  axis  orien¬ 
tation  is  constant  at  the  value  prescribed  in  (E-l)  for  a  line  through 
the  center  of  the  segment. 

The  segmentation  of  the  ring  has  two  consequences.  Equations  (E-2) 
and  (E-3)  are  multiplied  by  form  factors,  depending  on  the  exact  shape 
of  the  segments.  For  trapezoidal  segments,  the  form  factor  is 

Nn 

„  _ N  it  S1n  “A 

KN  "  cos  *-R —  (E-4> 

where  M  is  the  number  of  pieces  around  the  ring.  For  a  reasonably  large 
value  of  M,  kn  Is  close  to  I,  e.g.  for  N  ■  1,  a  dipole  magnet,  and  M  *  8 
pieces,  KN  »  0.90.  The  second  consequence  of  segmentation  is  that  the 
field  distribution  Inside  the  ring  Is  no  longer  a  perfect  multi  pole,  but 
contains  harmonic  field  errors.  The  fields  in  a  multipole  ring  may  be 
expressed  as 


B*(z)  -  Bx  -  1By  *  Z  bn(  ~  )n_1  (E-5) 

where  the  bn's  are  the  harmonic  multipole  field  coefficients.  It  can  be 
shown7  that  In  a  perfect  segmented  multipole  ring,  only  certain  terms 
appear  in  the  series  in  (E-5): 
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n  =  N  +  vH,  v  =  0,1,2,  ...  (E-6) 

For  a  dipole  with  M  =  8,  this  series  is  n  =  1,9,17,...  and  the  first 
error  harmonic  is  the  n  =  9,  18-pole.  The  strength  of  this  error  depends 
on  the  exact  geometry  of  the  pieces  and  the  ratio  (r1/r2)  .  For  8 
trapezoidal  pieces  in  a  dipole  with  (r1/r2)  =  ^  ,  | b 9/b x |  *  0.095. 

The  flux  distribution  in  such  a  ring  is  shown  in  Figure  8. 


2.  Helical  Wiggler:  An  REC  permanent  magnet  helical  wiggler  can 
be  constructed  by  an  axial  stack  of  dipole  rings  where  each  ring  is 
rotated  about  the  axis  relative  to  its  neighbors.  For  the  dipole  ring 
geometry  of  Figure  7  the  expression  for  the  amplitude  of  the  rotating 
dipole  field  on  axis  in  such  a  wiggler  1sa 


sin  y 

B„  *  Br  Cl— [T(X1)  -  T(x 2) ] 
J 


with  the  following  definitions: 

.  2* 
sin  ~M 

Ci  *  — ^ —  »  the  form  factor  for  circular  segments, 


(E-7) 


J  =»  the  number  of  slices  per  period  X,  and 


x  2ltri 

T(x)  «  K0(x)  +  j  Ki(x),  xi  » 


where  K0  and  are  modified  Bessel  functions  of  the  second  kind.  The 
function  T(x)  Is  given  In  Figure  9.  Both  Ci  and  the  factor  depending  on 
the  number  of  slices  per  period  will  be  close  to  1  for  reasonable  designs. 
T(x2)  can  be  made  small  for  reasonable  values  of  x2  »  2 wr2/X.  Therefore 
the  dominant  term  which  determines  the  helical  wiggler  strength  will  be 
the  value  of  T(x  1 )  corresponding  to  the  ratio  r^/X.  Although  (E-7)  is 
strictly  correct  only  for  the  case  of  circular  segments  in  the  dipole 
ring.  It  Is  a  good  approximation  to  the  solution  with  alternate  shaped 
segments  such  as  trapezoids  and  rectangles,  if  the  proper  form  factors, 

Ki,  are  substituted  for  CL. 
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COMPOSED 


TABLE  III 


r  2 

Segment 

Shape 

Bi  (Kgauss) 

b9 

F7<*> 

2 

trapezoid 

5.42 

9.5 

2 

rectangle 

4.18 

0.2 

3 

rectangle 

5.66 

0.01 

A  preliminary  analysis  of  a  helical  wiggler  design  with  = 

2  kilogauss  and  \  =  5  cm  has  been  made.  The  first  step  was  to  investi¬ 
gate  different  configurations  for  the  dipole  ring.  Table  3  gives  the 
results  for  infinitely  long  (2D)  rings  assuming  Br  *  8.7  kilogauss. 

For  mechanical  reasons,  a  design  with  rectangular  pieces  is  preferable 
as  this  allows  for  mechanical  support  of  the  pieces  and  the  possibility 
of  moving  the  pieces  to  tune  out  harmonic  errors,  as  will  be  discussed 
below. 

Figure  10  shows  a  helical  wlggler  design  concept  with  the  following 
parameters. 

sin  11 

-  (8.7)  (.9)  (.9) 

7 


xi  « 


2*  (1.5) 


1.85 


T(xi)  -  .3 

2*  (4.51 

X2  ■  "  1 g  '  *  5.65 

T(x2)  »  0.01 


This  gives  a  wlggler  amplitude  of 

Bh  -  (8.7)  (.9)  (.9)  (.29)  -  2.04  Kgauss 


3.  Dipole  Ring  Tuning:  Although  many  REC  multiple  magnets  have 
been  built  with  no  provisions  for  tuning  the  multipole  field  errors,  the 
requirements  on  the  uniformity  of  the  magnet  pieces  are  very  stringent 
for  high  precision  field  distributions.  A  way  to  understand  how  a 
tuning  procedure  may  be  implemented,  is  to  consider  the  field  produced 
by  an  individual  magnet  block.  In  general  a  single  block  produces  all 
harmonics  which  may  be  expressed  as 

B*  -  iBj  =  Z  C  z"'1  (E-8) 

y  n-1 

Let  us  call  the  block  bisected  by  the  x-axis  in  a  multi  pole  assembly  the 
reference  block,  as  shown  in  Figure  11.  The  Cn  in  (E-8)  are  produced  by 
this  block.  The  resultant  fields  produced  by  an  assembly  of  M  blocks  in 
a  symmetrical  array  are  just  the  sum  of  the  harmonics  produced  by  a 
single  block,  i.e.  bn  in  (E-5)  Is  equal  to  MCn.  Of  course  only  certain 
n's  are  allowed  according  to  (E-6)  all  other  harmoni.cs  exactly  canceling. 
The  Cn  for  an  arbitrarily  shaped  reference  block  may  be  derived  by 
evaluating  the  following  integral 

c  ,  Ik.  Ji z*  (E-9) 

4 it  i  zn 

where  Br  is  the  complex  representation  of  the  residual  induction.  If 
the  position  of  the  reference  block  is  perturbed,  the  change  in  the  Cn 
can  be  calculated  directly  from  (E-9).  The  result  is 

ACn  ■  -n  Az  Cn+1  (E-10) 

The  same  relationship  exists  for  every  block  in  the  assembly  with  a 
different  absolute  phase  angle  for  each  block. 

The  consequences  of  the  analysis  above  are  that  by  small  pertur¬ 
bations  of  the  position  of  the  REC  blocks  In  a  multipole  assembly, 
specific  harmonic  errors  of  given  amplitude  and  phase  (real  and  Imaginary 
part)  may  be  produced. 

\ 

A  particular  multi  pole  tuning  procedure  then  could  be  implemented  as 
follows.  First,  measure  the  amplitude  and  phase  of  each  harmonic 
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multipole  field  error.  This  can  be  done  with  harmonic  multipole 
measuring  equipment  as  described  in  reference  9.  Second,  calculate  the 
set  of  magnet  block  radial  perturbations  which  create  the  negative  of 
the  measured  errors  with  the  proper  phase.  Depending  on  the  criterion 
which  is  chosen  to  reduce  the  field  errors,  one  may  wish  to  eliminate 
exactly  a  number  of  low  order  errors,  or  the  strength  of  a  larger  number 
of  multi  pole  errors  may  be  reduced  by,  for  example  a  weighted  least 
squares  minimization  procedure.  In  general  it  requires  at  least  the 
radial  movement  of  two  magnet  blocks  to  cancel  one  harmonic  error,  since 
the  bn  are  complex. 

F.  Accelerator  Studies 

The  success  of  the  two-stage  FEL  experiment  will  depend  on  the 
proper  functioning  of  all  equipment  in  the  system.  This  includes  the 
accelerator  that  supplies  the  beam.  Beam  current,  pulse  length,  and 
voltage  stability  are  critical  factors  in  the  experiment  that  are  deter¬ 
mined  by  the  accelerator,  and  the  amount  of  power  that  can  be  obtained 
from  the  FEL  is  also  limited  by  the  accelerator. 

During  the  period  covered  by  this  report  we  have  studied  the 
adequacy  of  the  present  UCS8  accelerator  system  for  the  two-stage  FEL 
experiment.  As  a  result  of  this  work,  we  have  identified  areas  in  which 
accelerator  improvements  are  necessary  to  Insure  the  success  of  the  two- 
stage  FEL  experiment  and  the  usefulness  of  the  experimental  results  for 
extrapolation  to  higher  laser  powers.  We  have  also  developed  a 
conceptual  design  for  a  very  high  power  FEL  and  identified  technologies 
needed  for  this  high  power  design  that  could  be  tested  using  the  UCSB 
accelerator. 

1.  Deficiencies  of  the  Present  UCSB  Accelerator;  In  the  present 
UCSB  electrostatic  accelerator, 10  the  charge  on  the  dome  of  the  accel¬ 
erator  Is  provided  by  a  pelletron  charging  chain  capable  of  supplying 
500  pamp  of  Input  current.  Initially  the  accelerator  will  operate  with 
a  2-amp  pulse.  Electrons  leaving  the  accelerator  pass  through  the  FEL 
amplifier  and  are  directed  back  up  the  accelerator  column  where  they  are 
recollected  and  reused.  The  charging  current,  therefore,  Is  only  needed 
to  start  up  the  accelerator  and  to  compensate  for  system  losses  during 
operation. 
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If  99.975%  of  the  electrons  leaving  the  dome  of  the  accelerator 
were  recollected,  an  average  current  of  2  amp  could  be  achieved.  It  is 
not  expected,  however,  that  collection  efficiency  will  be  this  good. 
Therefore,  the  accelerator  will  be  operated  in  a  pulsed  mode.  As  a 
pulse  of  constant  current  is  emitted  from  the  accelerator,  the  charge  on 
the  dome  will  be  depleted  and  the  accelerating  voltage  will  begin  to 
fall  or  droop.  The  length  of  the  pulse  will  be  a  function  of  the  droop 
that  can  be  tolerated  in  the  FEL.  Droop  voltage  for  a  2-amp  beam  in  the 
UCSB  accelerator  is  plotted  as  a  function  of  pulse  length  for  a  number 
of  assumed  current  loss  rates  in  Figure  12.  For  a  single-stage  long- 
wavelength  FEL  several  kilovolts  of  droop  could  be  tolerated  before  the 
laser  ceases  to  produce  gain  at  the  desired  frequency.  For  the  second 
stage  of  a  two-stage  FEL,  the  maximum  droop  that  could  be  tolerated 
would  be  a  few  hundred  volts.  Even  if  a  voltage  droop  could  be 
''tolerated",  interpretation  of  the  data  from  an  experiment  with  voltage 
droop  would  be  very  complicated  and  difficult  to  extrapolate  to  high 
power  operation. 

The  length  of  the  electron  pulse  Is  determined  by  the  rate  at  which 
charge  Is  lost  by  the  dome  of  the  accelerator.  For  the  present  system  a 
pulse  length  of  75  usee  ,  determined  by  the  amount  of  droop  that  could 
be  tolerated  in  the  first  stage,  has  been  predicted  for  operation  at  2 
amp.  If  the  electron  current  Is  Increased  by  a  factor  of  10,  it  is 
anticipated  that  the  pulse  length  will  be  decreased  by  at  least  this 
factor. 

In  a  free  electron  laser  oscillator,  a  finite  amount  of  time  is 
required  for  the  laser  pulse  to  build  up  ivom  spontaneous  emission  to  a 
saturation  value  at  which  the  gain  in  the  laser  just  covers  optical 
system  losses.  The  time  required  for  this  buildup  in  the  Stanford  FEL 
Is  a  minimum  of  a  few  tens  of  microseconds.  For  both  the  second  stage 
laser  field  and  the  first  stage  pump  field  to  build  up  from  spontaneous 
emission  will  probably  require  at  least  several  tens  of  microseconds. 

It  Is,  therefore,  not  certain  that  with  the  present  accelerator  design 
the  electron  pulse  will  be  long  enough  to  permit  the  laser  pulse  to 
reach  Its  predicted  saturation  value. 

We  assume  that  a  maximum  effort  will  be  made  to  lengthen  the 
electron  pulses  by  reducing  beam  losses.  The  predicted  75  usee  pulse 
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Figure  12 


length  for  a  2-amp  beam  and  7.5  nsec  pulse  length  for  a  20-amp  beam  are 
obtained  by  assuming  these  efforts  will  be  successful.  In  order  to 
increase  the  pulse  length  further,  it  will  be  necessary  to  correct  for 
the  voltage  droop  that  limits  the  pulse  length. 

2.  Control  of  Voltage  Droop:  The  rates  for  charging  and 
discharging  the  dome  of  the  accelerator  are  given  by 

lr  I, 

Kc  a  r"  and  K0  *  7T  (F_1) 

where  is  the  charging  current,  1^  is  the  current  that  is  not  recol¬ 
lected  by  the  accelerator,  and  C  is  the  capacitance  of  the  dome.  The 
net  rate  of  discharging  is  given  by 


where  AV  is  the  voltage  drop  of  the  terminal.  For  the  UCSB  accelerator 
current  losses  during  the  electron  pulse  are  much  greater  than  the 
charging  current,  i.e-.  IL  »  1^. 

Without  correction  for  voltage  droop,  the  electron  pulse  length  is 
determined  by  the  voltage  droop  that  can  be  accepted.  This  discharge 
time  will  be  given  by 
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The  time  required  to  recharge  the  terminal  between  pulses  is  given  by 


At 


C 


( F— 4) 


and  the  duty  cycle  is 


D 


At. 


+  At. 


(F-5) 


For  example,  if  Iq  ■  500  pamp  ,  IL  »  20  mamp  and  C  3  200  pF,  then 
Kc  *  2.5  V/psec  and  Kg  ■  100  V/psec  .  If  we  allow  a  AV  of  10  kV  then 
the  pulse  length  A£q  is  100  usee  ,  the  recharging  time  is  4000  psec  and 
the  duty  cycle  is  0.025.  This  is  shown  schematically  in  Figure  13(a). 
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The  voltage  droop  can  be  eliminated  by  providing  a  voltage  ramp  to 
accelerate  the  electrons  by  exactly  the  amount  that  the  terminal  voltage 
is  reduced  (Figure  13b).  This  voltage  ramp  would  be  provided  by  a  power 
supply  in  the  dome  of  the  accelerator  that  ramps  the  voltage  of  the 
electron  gun  and  electron  collector  relative  to  the  dome  during  the 
current  pulse.  The  power  needed  to  operate  the  power  supply  would  be 
provided  by  the  generator  in  the  dome.  Any  change  in  beam  focusing  pro¬ 
duced  by  the  accelerating  gap  could  be  corrected  using  variable  voltage 
Einzel  lenses. 

When  the  voltage  droop  has  been  cancelled,  the  pulse  length  is  no 
longer  determined  by  the  maximum  droop  the  FEL  can  accept,  but  by  the 
maximum  droop  the  ramped  power  supply  can  correct.  The  latter  may  be 
greater  by  a  factor  of  a  few  than  the  former.  The  size  of  the  power 
supply  and  the  breakdown  voltage  for  the  accelerating  gap  limit  the 
droop  that  can  be  corrected  to  about  30  keV.  If  the  charging  current 
remains  constant,  there  will  be  no  change  in  duty  cycle,  so  that  the 
time  required  to  recharge  the  terminal  increases  as  the  electron  pulse 
length  is  increased  (Figure  13c).  Increasing  the  charging  current  will 
have  practically  no  effect  on  the  electron  pulse  length  when  IL  »  Iq, 
but  will  increase  the  duty  cycle  (Figure  13d). 

3.  Control  of  Voltage  Ripple:  The  pelletron  chain,  which  Is  used 
to  transfer  charge  to  the  HV  terminal,  is  charged  by  induction.  It  can 
be  charged  negatively  at  the  accelerator  base,  and  positively  at  the 
terminal  to  double  the  net  charge  transferred.  Ripple  is  caused  by 
imbedded  charge  which  is  not  transferred  in  the  charge-discharge  process. 
The  amount  of  imbedded  charge  can  vary  with  time.  The  ripple  voltage 
has  two  main  frequency  components:  one  at  400  Hz  arising  from  the 
charging  variations  of  Individual  pelletron  links,  and  one  at  a  much 
lower  frequency,  ~  3  Hz,  corresponding  to  the  complete  chain  cycle 
time.  For  a  completely  unregulated  machine,  the  amplitn^0  of  the  ripple 
Is  on  the  order  of  kilovolts. 

Conventional  stabilization  methods  utilize  a  coarse  and  fine  control 
to  minimize  ripple.  The  coarse  control  can  only  control  ripple  to  about 
1  part  in  103  of  the  total  accelerator  voltage,  whereas  the  fine  control 
usually  reduces  the  ripple  to  1  part  in  10^  or  less.  The  state-of-the 
art  in  ripple  control  is  a  few  parts  in  105,  Coarse  control  for  the 
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UCS8  machine  utilizes  a  generating  voltmeter  to  detect  voltage  changes 
on  the  terminal  and  a  corona  discharge  to  modulate  a  leakage  charge  for 
correction. 

For  fine  control  of  ripple  the  signal  from  an  energy  sensitive 
diagnostic,  such  as  a  stabilized  analyzing  magnet  and  slit  detector, 
could  be  used  to  modulate  the  voltage  on  the  same  power  supply  already 
being  used  to  correct  for  voltage  droop. 

4.  Accelerator  Power  System  Limitations:  Although  the  20  kw 
generator  and  associated  motor  and  power  supplies  in  the  (JCSB  accel¬ 
erator  are  adequate  at  present,  they  may  not  be  able  to  provide  adequate 
power  when  the  beam  current  is  increased.  If  1%  of  the  electron  energy 
is  converted  to  photon  energy  in  a  2-amp,  3-MeV  beam,  the  instantaneous 
power  that  must  be  supplied  is  60  kw.  For  a  20-amp  beam  this  will 
increase  to  600  kw.  Because  the  gain  in  both  the  pump  and  laser  fields 
will  be  higher  in  the  20  amp  case,  more  than  1%  of  the  electron  energy 
may  be  converted  to  photon  energy  when  the  system  reaches  its  equilibrium 
operating  point.  In  order  to  accommodate  this  higher  energy  loss  it 
will  be  necessary  to  Increase  the  voltage  on  the  electron  gun  and 
Increase  the  voltage  spread  that  the  collector  can  accept.  Correction 
for  voltage  droop  will  permit  increasing  the  electron  pulse  length  and 
hence  the  time  over  which  peak  power  is  drawn  from  the  power  supplies  in 
the  dome.  This  will  Impact  the  amount  of  energy  that  must  be  stored  to 
provide  for  peak  power  loading. 

5.  Design  of  an  Accelerator  for  a  High  Power  FEL:  Even  if  the 
voltage  stability  of  the  UCSB  accelerator  is  improved  it  will  not  be 
capable  of  driving  a  high-average-power  FEL.  One  reason  for  this  Is 
that  the  total  power  that  can  be  generated  In  both  the  pump  and  the 
laser  fields  of  the  FEL  cannot  exceed  the  generator  power  In  the  dome  of 
the  accelerator.  In  the  UCSB  device,  there  Is  only  a  20  kw  generator  in 
the  dome  of  the  accelerator. 

Other  factors  also  limit  the  average  power  that  can  be  obtained 
using  the  UCS8  accelerator  design.  When  operating  in  a  pulsed  mode,  the 
laser  pulse  must  build  up  from  spontaneous  emission  for  each  pulse. 
Therefore,  the  laser  will  at  best  operate  at  Its  saturation  power  value 
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for  only  a  fraction  of  the  time  the  electron  pulse  is  on,  and  the  elec¬ 
tron  pulse  will  be  on  only  a  small  fraction  of  the  time  between  pulses. 
The  low  duty  cycle  is  the  result  both  of  inadequate  charging  current  to 
permit  continuous  operation  and  inadequate  power  to  provide  energy  to 
the  electrons  for  continuous  operation. 

At  KMSF  we  have  developed  a  design  concept  for  an  electrostatic 
accelerator  that  could  provide  both  the  input  current  and  the  power 
needed  to  operate  a  high-average-power  FEL.  In  the  accelerator  design 
we  are  considering,  an  insulated  core  transformer  (ICT)  is  used  to 
provide  the  input  current  to  the  dome  of  the  accelerator.  An  ICT  could 
provide  a  few  tens  of  milliamperes  of  input  current  at  voltages  of  a  few 
MeV  directly  to  the  dome  of  the  accelerator  with  90%  wallplug  efficiency. 
A  diagram  of  an  accelerator  incorporating  an  ICT  is  shown  in  Figure  14. 

To  provide  adequate  power  for  operating  the  FEL,  four  samarium 
cobalt  permanent  magnet  generators,  each  of  which  is  rated  at  250-400  kw 
and  weighs  120  lb.,  are  located  In  t.ie  dome  of  the  accelerator.  The 
generators  are  powered  by  two  800  hp  motors  with  counter-rotating  drive 
shafts  located  outside  the  accelerator.  The  Insulating  drive  shafts  may 
be  operated  directly  using  a  ferro-fluidlc  feedthrough,  as  shown  In 
Figure  14,  or  using  a  fluid-driven  turbine  to  transmit  power  to  the 
drive  shaft  inside  the  accelerator.  The  permanent  magnet  generators 
operate  at  20,000  rpm  and  a  gear  box  would  be  used  to  operate  the  drive 
shaft  at  2000  rpm.  The  overall  efficiency  of  producing  power  in  the 
dome  of  the  accelerator  is  estimated  to  be  about  80%. 

All  of  the  equipment  needed  to  build  this  accelerator  Is  commer¬ 
cially  available.  The  major  technological  problem  that  needs  to  be 
demonstrated  is  efficient  recollection  of  the  high  current  electron 
beam.  Also,  operation  of  the  generators  and  the  entire  power  transmis¬ 
sion  assembly  must  be  demonstrated  Inside  an  electrostatic  accelerator. 

The  voltage  of  the  ICT  has  a  periodic  ripple  of  about  1%  at  360  Hz. 
This  ripple  could  be  reduced  by  a  factor  of  1000  using  the  filter  circuit 
shown  In  Figure  15.  It  Is  necessary  In  constructing  this  filter  that 
the  accelerator  and  ICT  column  domes  remain  separate.  This  is  shown  in 
Figure  14.  Further  correction  for  small,  nonperiodic  voltage  ripple 
could  be  accomplished  with  a  beam  energy  analyzer  and  feedback  system. 
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Ferro-Fluidi 
Feedthrough  ( 


riWiTA% 


ICT  FILTER  CIRCUIT 


ASSUMPTIONS 

*  360  Hz  (3-Phase,  FWR)  with  1%  Ripple 

*  60  Hz  with  0.02%' Ripple 

*  Cj_  =  C2  =  200 /up.? , 

*  RL  **  3500  mQ  (1  mA  bleeder  current) 

RESULTS 

*  ^FILTER  *  200  MQ 

*  Ripple  (360  Hz)  is  1  in  105 

*  Ripple  (60  Hz)  is  1  in  105 

*  Power  Dissipation  is  30  mA  x  35  kV  *  1  kW 


Figure  15 


As  an  example,  we  consider  an  accelerator  with  3.5-MeV,  20-amp  DC 
beam  that  is  recirculated  with  99.9%  charge  collection  efficiency. 

The  ICT  makes  up  the  20-mamp  current  loss.  If  an  average  of  1.4%  of 
the  electron  beam  energy  is  converted  to  photon  energy  in  the  FEL, 

1  MW  of  power  must  be  supplied  to  the  recollected  electron  beam  by  the 
generators  in  the  dome  of  the  accelerator.  Since  the  ICT  is  only  90% 
efficient  and  the  generator  system  is  ~  80%  efficient,  heat  will  be 
generated  in  the  accelerator  tank  that  would  have  to  be  removed  by  an 
active  cooling  system. 

Although  the  accelerator  column  and  ICT  are  shown  in  tandem  in 
Figure  14,  it  would  also  be  possible  to  locate  them  along  side  each 
other  for  a  more  compact  design  (Figure  16). 

G.  Multidimensional  FEL  Simulation  Code 

The  one  dimensional  resonant  particle  and  multiparticle  simulation 
codes  that  have  been  developed  at  KMSF  have  enabled  us  to  obtain  general 
scaling  relations  for  the  two-stage  FEL  as  well  as  more  detailed  infor¬ 
mation  about  optimization  of  laser  gain  and  conversion  efficiency  during 
the  pulse  buildup.  These  codes  cannot,  however,  provide  information 
about  the  mode  structure  of  the  optical  beam,  beam  quality,  or  output 
coupling,  or  provide  Information  about  the  actual  electron  trajectories 
In  the  Interaction  region.  The  1-0  codes  do  not  take  the  radial  profile 
of  the  electron  beam  and  optical  beam  into  account,  and  therefore  provide 
only  approximate  values  for  laser  gain  and  conversion  efficiency,  which 
are  not  sufficiently  accurate  for  detailed  comparison  with  experiment. 

To  Improve  our  capability  to  model  a  two-stage  FEL  experiment  we 
have  developed  two  and  three  dimensional  computer  simulation  codes  and  a 
two  dimensional  resonant  cavity  design  code.  In  this  section  we  will 
present  the  general  three  dimensional  formalism  for  propagation  of  the 
laser  field  and  the  electromagnetic  pump  field  in  the  FEL.  The  method 
of  solution  of  the  wave_ equation  using  fast  Fourier  transforms  In  a 
Cartesian  coordinate  system  will  be  described,  and  then  It  will  be  shown 
how  this  same  formalism  can  be  applied  to  a  system  in  which  cylindrical 
coordinates  are  employed. 


In  the  development  of  the  multidimensional  codes  a  number  of  subtle 
points  had  to  be  considered  to  insure  that  the  codes  would  be  both  fast 
and  accurate.  In  particular,  it  had  to  be  shown  that  the  incremental 
solution  of  the  wave  equation  would  be  unitary,  so  that  no  artificial 
gains  or  losses  would  be  introduced  in  the  absence  of  a  driving  source 
term.  To  solve  the  wave  equation  in  cylindrical  coordinates  using  fast 
Fourier  transforms,  a  transformation  is  made  on  the  radial  part  of  the 
wave  equation. 

Since  the  experiment  we  are  planning  will  have  cylindrical  symmetry, 
a  2-0  code  has  been  written  which  assumes  cylindrical  symmetry.  To  use 
the  code  for  design  optimization  of  an  FEL  amplifier  it  is  necessary 
to  know  what  effect  varying  an  axial  electric  field  (or  equivalently 
tapering  the  period  of  a  wiggler  magnet)  would  have  on  the  resonant 
phase  of  the  ponderomotlve  wave  as  a  function  of  radial  position.  This 
is  determined  by  developing  a  description  of  the  electron  dynamics  for 
the  resonant  particle  In  the  2-0  code. 

1.  The  Three  Dimensional  Wave  Equation;  The  propagation  of  an  electro¬ 
magnetic  wave  in  space  Is  determined  by  the  three  dimensional  Maxwell 
wave  equation 

'  •'*«'’  (s'1) 

We  will  solve  this  equation  for  the  general  case  in  which  both  a 
laser  field  Aj_  and  an  electromagnetic  pump  field  Ap  combine  to  form  a 
ponderomotlve  force  which  acts  on  the  electrons. 

The  currents  used  In  the  solution  of  the  wave  equation  In  any  time 
Interval,  At,  are  obtained  from  the  motion  of  the  electrons  during  that 
time  Interval.  The  solution  of  the  wave  equations  In  the  1-D  case  was 
made  simpler  by  using  two  assumptions:  1)  the  slowly  varying  field 
approximation  which  allowed  linearization  of  the  left  hand  side  (LHS)  of 
the  wave  equation,  and  2)  the  proper  choice  of  propagation  vectors  on 
which  to  project  the  wave  equation.  The  latter  allowed  the  separation 
of  the  differential  equations  for  the  amplitude  and  phase  of  the  fields 
and  simultaneously  provided  exact  analytic  forms  for  the  single-particle 
currents  (the  RHS  of  the  wave  equation).  It  will  be  shown  that  using 


these  assumptions  the  wave  equation  can  be  written  in  the  form 


[7i2  +  2ikL  (-gj  +  AL  =  -M-oJ 


(G-2) 


where  ki  is  the  longitudinal  component  of  the  wave  number  for  the  laser 
field  and  J j_  is  the  perpendicular  component  of  the  electron  current. 
Likewise  for  the  pump  field 


[vi2  +  2ikp  (<g|  -  3-g|)]  Ap  =  -no«Ji 


(G-3) 


In  the  3-0  case,  we  generalize  the  vector  potentials  to  have  ampli¬ 
tude  and  phases  that  depend  on  x,  y,  z  and  t  but  which  propagate  only  in 
the  z-dl recti  on.  The  latter  property  suggests  we  find  a  way  to  represent 
the  potentials  used  in  the  wave  equation  so  that  the  evolution  of  the 
amplitude  and  phase  part  can  be  handled  separately  from  the  propagation 
part. 

The  vector  potentials  for  the  laser  field  and  the  electromagnetic 
pump  field  are  given  by 

Al  3  AL(x,t)  [x  cos(kLz  -  aj|_t  +  $L)  -  y  sin (kLz  -  i^t  +  $|_)]  (G-4) 

and 

+  ♦  *  - 

Ap  ■  -Ap(x,t)  [x  cos(kpz  +  <Dpt  +  $p)  +  y  sin(kpz  +  ojpt  +  4>p)]  (G-5 ) 

where  \  3  \(x,t),  $p  *  $p(x,t). 

We  now  rewrite  these  in  a  form  suitable  for  separating  the  z,t  propa- 


Note  that 


i  (k.  Z-U>  t)  -  r  *  A  _ 

-  (x,iy)  =  Lx  cos(kLz-o)Lt)  -  y  sin (k^z-o^t )  J 

+  i  [+  x  sin(kl_z-(JLt)  +  y  cos(kLz-W|_t)  ] 


=  e i  -  ie? 
XL  ZL 


where 


A  A 


and 


Also  note  that 


where 


e^  s  x  cos(kLz-Gjj_t)  -  y  sinCk^z-c^t) 


e2L  2  "x  s1n(kLz-ojj_t)  -  y  cos(kLz-o)Lt) . 


i(kpz+upt) 


(x,-iy)  •  -(elp  -  1 e2p ) 


A  A 


and 


elp  5  “x  cos(kpZ+upt)  -  y  s1n(kpz+upt) 


e2p  5  x  sin(kpz+wpt)  -  y  cos(kpz+wpt) , 


1*i 


We  further  define 

Aj_  =  AL(x,t)e 

Ap  5  Ap(x,t)e1,>P 


*  A 


fL  2  elL  -  ie2L 


(G— 9 ) 


(G— 10) 


( G— 11) 


(6-12) 


(G— 13) 


(G-14) 


(G-15,16) 


fp  s  -(elp  -  ie2p)  (6-17,18) 


Aj  »  A>  carry  the  amplitude  and  phase  shift  Information  which  depend  on 

'  a  A 

(x,t)  and  the  vectors  f ^ •  fp  depend  only  on  (z,t).  Then 


♦  A  -  A 

Aj_  ■  Re [£^_f ^ +  c.c. 


(6-19) 


and 
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♦  r  A 

Ap  -  -Re 1 Ap  f p 1  -  - 

V2Apfp  +  c*c. 

(G-20) 

Defining 

“L  =  (ktz  - 

“Lt), 

(G-21) 

and 

*o  s  \  +  ®p  3  (kL+kp)z  - 

(G-22) 

it  can  be  shown  that 

"el£  elp  3  e2L*e2p  3  cos 

®1l  ®il  =  ll  etc- 

(G-23,24) 

el|[  e2p  3  e2£  elp  3  s"*n  7° 

A  A  A  A 

ei  •  e?  =  ei »  eo  =  0. 
iL  dL  lp  zp 

(G-25,26) 

It  can  also  be  shown  that 

A  A  * 

Vfl  >2, 

A  A  + 

VfP  =2 

(G-27,28) 

A  A  X 

A  *  * 

VfP  -0, 

Vfl  ■» 

(G-29,30) 

A  ★  A  * 

-ilo 

and  fp  *f^  *  2e  , 

(G-31) 

Also, 

a  2  “2  A  *  2 

A  ★  2 

fL  3  fP  3  <f|_  ) 

3  (fp  )  3  o. 

(G-32) 

The  result  of  this  formalism  is  that  the  complex  form  of  the  fields  can  be 
written  as  projections  of  the  real  vector  potentials  onto  the  complex 
vectors : 

•+  »  *  *►  *  * 

»  A^»f^  ,  jip  »  -Ap»fp  .  (G-33,34) 


Slnce-sf  <*‘>L,P  '  kl,P  <*2>L,P  and  at  <«2>U,P 


-k 


UP 


(e  i) 


L,P 


then: 
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v 


(G-35,36) 


Simi 1 arly 


T5F  fl  -  ikLfL  • 


bz  fP  =  lkPfP 


3  f 
■St  fL 


_a£ ;  -  ir  2  f 

az2  fL  -  -kL  ' L 

_a£  f  -  _k2  f 

bz 2  TP  kP  P 


’HV  ^2 


■St  fP  =  iuJPfP  *  ^72 


*2  *  a 
a?  fL  "  _UJL  fL 
a2  ;  2  ; 

at2  fP  =  _wp  fP* 


(G-37,38) 


(G-39,40) 

(G-41,42) 


Evaluation  of  the  LHS  of  the  Wave  Equation 

Next  we  evaluateQ2A  for  each  field,  where 

□  2  =  [7  S+_Si.ij_4]. 

1  az2  c2  at2  J 


For  the  laser  field 

a2  l  a2i  e\ 


,  .  a2A.  aA.  af,  ,  a2f, 

r  o  i  Q"i  /-l  .  *  i  r  1  2  -L  .  -L  L  .  1  .  _ Li 

^bP  "  C7  at^  ^2  4fi^  fL  aP”  +  az  az  +  2^.~bP~* 


(G— 43) 


,  ,  •  a2A,  aA.  af,  .  a2f,  , 

'  c7  $  f  L~aF  +  %  “at  +  \  4“at^ 

.  8A.  .bA 

*  1kifL  t-i? +  i  -S’]  (G-44) 

a2A.  a2A, 

where  we  have  assumed  k.  =  «,  /c  and  — -jp-  -  — -tr  "  °»  that  is, 

u  n.'  az-:  at^ 

the  fields  are  slowly  varying  in  the  propagation  direction.  Then 

a2\  -  [i  7^  *  1kL  (jf  *  ±  $]  4?l  *  c.c.  -  -MO,  (6-4S) 


or 


R«  {[{  7^  ♦  (kL(^  ♦  Z  3f)]  VL)  *  -  !*  Jr  (S'46> 
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4 


°i 5 1?  fi5  *  «L  (*£♦?*>] 


ala* 


(6-47) 


then 


(DL4)fL  +  (D^)  fL  = 

To  eliminate  from  the  LHS,  project  both  sides  onto  4*: 

(0l4.)fW  +  (DL4}*(fV)2  3 


(G-48) 


(G-49) 


or 


Di4  = 


Similarly,  for  the  pump  field: 


-1 


Re  (Cf  V*  a/,) -fj,. 


and 


Dpip  *  f1  VfP*  • 


(G-50) 

(fi-51) 

(G— 52) 


Evaluation  of  the  RHS  of  the  Wave  Equation 
We  look  now  at  the  form  of  the  currents 

»  f  ecpx.  (G-53) 

where  + 

Pll*  (  m^  “  m£  )  »  (G~54) 

A  *  Al  +  Ap,  (G-55) 

and  the  sum  Is  taken  over  all  the  electrons  in  a  ponderomotl ve  wavelength. 
Using  (G-55),  (G-33),  (G-20),  (G-29),  and  (G-31)  It  can  be  shown  that 

*  -  *  *  -1*0 

A  •  fL  -  4  -  Ap  e  (G-56) 
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Simi larly 


A  •  fD  =  A,  e 


*  -1*o 


~  Ap  • 


(G-57) 


With  these  results 


and 


Ji‘fL 


VfP 


a  r*  /'  I  i  *  ★  p  f  ★  *1  «  . 

•i  I  c— 77  •  fL  -  77  [4  “  Ap  e  ] ) 


■  |  ^  •  V  *  y7  [-Ap  +  4*e  1701  ])• 


(G— 58) 


(G-59 ) 


Finally,  the  linearized  wave  equations  for  the  pump  and  laser  fields  are 
given  by 

*  -iYOf 


[£  V  +  1kl^5z  +  c  at)  1  4*  ‘  /  m  E  C“7T  *  fL*  -  7?  C^Li  ' 


)] 


and 


7j_2  +  ikp(-g|  -  i  af)  ]  Ap  -  5  l  •  V  -  Y?  ^i  +  ^  1101  5] 

2.  Incremental  solution  of  the  wave  equation 

By  transforming  to  a  coordinate  system  moving  with  the  wave,  the 
time  derivation  in  the  linearized  wave  equation  (G-2)  can  be  eliminated. 
The  linearized  wave  equation  can  then  be  written  in  the  general  form 

9A 


OA  1  w  2  a  c 

TSz  *  IF  7i  A  -  S 


(G— 62) 


where  S  is  a  source  term.  A  separate  equation  is  needed  for  each 
transverse  component.  An  incremental  solution  to  (G-62)  is 

A (z+dz)  ■  A(z)  +  7J.2Mz)  ■  AzS 

■  (l  +  ?l2)  a(z)  -  AzS 


(6-63) 


Equation  (G-63)  is  not  unitary.  An  incremental  solution  to  (G-62)  which 
is  unitary  is 


1**  7  2 

A(z+Az)  *  e  2k  1  A(z)  -  AZS 


(G-64) 


(G— 60 ) 

(G— 61 ) 
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Equation  (G-64)  can  be  solved  in  Cartesian  coordinates  using  Fourier 
transforms.  If  F  is  the  Fourier  transform  operator  and  F"1  is  its 
inverse,  then 


7 

A ( z+Az )  =  e  **  1  F_1FA(z)  -  AzS 

2+k  2\ 

=  F"1  [e^  x  y  F  A (z)  ]  - 


(6-65) 


Using  a  computer  code  with  fast  Fourier  transform  subroutines, 

(6-65)  can  be  used  to  numerically  propagate  an  electromagnetic  wave. 

j  42  _  2 
“5F  7i 

The  source  term  S  drives  the  wave  and  the  operator  e  spreads  the 

3  V  2 

wave  due  to  diffraction.  To  understand  why  the  operator  e  2  1  is  used 

rather  than  the  operator  (1  +  v^2)  of  (G-63),  we  consider  the  case  in 

which  the  source  term  S  a  0.  In  this  case  there  can  be  no  gain  and  the 
Fourier  components  A(kx,ky,z)  =  F  A(x,y,z)  should  not  vary  as  a  function 
of  z.  This  property  holds  for  the  exponential  operator  since 


“2iT  (kx2+lcy2) 

|e  x  y  |  =  1 

It  is  not  true,  however,  for  the  operator  of  (G-63).  If 

A(kx,ky ,z+Az)  =  (l  +  (kx2  +  ky2))  A(kx,ky,z) 


(6-66) 


(6-67) 


I  A(kx  ,ky ,z+Az )  |  2  *  [l  +  (#  (^x2+ky2))2  ]  |A(kx,ky  ,z)  |  2  (6-68) 

Equation  (G-68)  says  that  each  mode  of  the  wave  will  be  amplified  by  a 

factor  1  +  jg-  (kx2+ky2)  on  each  iteration  of  the  computer  simulation. 

After  many  Iterations  the  wave  function  would  be  dominated  by  high 
frequency  noise,  since  the  largest  values  of  kx  and  ky  produce  the 
greatest  gain. 


.  solution  of  the  wave  equation  in  cylindrical  coordinates 

The  two-stage  FEL  system  we  are  studying  has  cylindrical  symmetry  and 
cylindrical  boundary  conditions.  It  would,  therefore,  be  convenient  to 
model  this  system  in  cylindrical  coordinates.  If  we  expand  A(x)  and 


S(x)  in  terms  of  their  angular  Fourier  components 

A(x)  =  l  a  .(r,z)  eU9 


(G-69) 


S(x)  =  I  s.(r,z)  e1 


JJ,3_OD 


and  substitute  these  into  (G-62)  we  obtain  a  series  of  equations 

a  -i  r  l  a  a  a2, 

3z  aA  2k  *■  r  3r  r  3r  +  "r^-'  aA  SA 


(G-70) 


( G— 71 ) 


The  symmetry  of  the  system  determines  which  values  of  A  are  needed  to 
describe  the  radiation  fields.  For  example,  for  a  cylindrically  symmetric 
system,  such  as  a  beam  with  ■»  Gaussian  Intensity  profile  in  an  open 
resonator  cavity,  only  the  A=0  equation  would  be  needed.  For  a 
TEoi  mode  in  a  cylindrical  waveguide,  only  the  Aai  equation  would  be 
needed  to  describe  the  wave  propagation.  We  cannot  use  the  fast  Fourier 
transform  (FFT)  technique  with  the  radial  part  of  the  Laplaclan  to  solve 
the  Incremental  equation  as  in  (G-65).  However,  it  is  possible  by  making 
a  change  of  variables  to  reformulate  the  wave  equation  in  a  way  that 
permits  the  use  of  fast  Fourier  transform  techniques  with  cylindrical 
coordinates.  If  we  define 


gx(r,z)  *  /ra*(r,z) 


(G-72) 


we  find  that  g  satisfies 


9*  *  2k  ^  ar2 


a2  a2-  V4 


]  gx  -  /r  Sj 


The  Incremental  solution  to  this  equation  is  given  by 

2k  W*  r2  * 

gA(r,z+Az)  *  e  9i(f,z)  -  Az/r  sA(r,z) 


(6-73) 


(G— 74) 


The  eigenfunctions  of  the  operator  (- T72  +  (A2-V4)/r2)  subject  to 


>•  --  i,':,  .7^1 


the  boundary  condition  gjj.(rmax)  =  0  are  given  by  /r  Jj^kp^r)  where  the 
JA's  are  ordinary  Bessel  functions.  The  knA  are  determined  by  the  boundary 
conditions  at  r^*,  that  is  J*(knJirmax)  =  0*  More  precisely 

[-  Is'75) 


In  order  to  utilize  the  FFT  method  we  rewrite  (G-74)  in  the  form 

i  AZ  (A2-  Vd  )  i Az  a2 
2k  Fr^  2k  ar7 

gA(r,z+Az)  -  e  e  g A( r ,z )  -  Az/r  s* 


-i  Az 
IT 
=  e 


-i  Az 

“21 r 


P2 


g^(r,z)  -  Az/r  s*  (G-76) 


where  the  operators  P2  and  V 

P2  =  - 


are  defined  by 
a2 


and 


V  = 


V4 

r2 


(G-77,78) 


The  error  resulting  from  assuming  the  operators  in  (G-76)  are  commutative 
Is  of  order  Az2,  which  can  be  made  arbitrarily  small. 


To  implement  (G-76)  we  must  restrict  the  radial  coordinate  to  a 
discrete  set  of  values  rj-jAr,  where  j*0,l,...,N-l  and  Ar*rmax/N. 

We  must  also  sample  gA(r,z)  at  these  same  points.  The  operators  in  the 
exponentials  must  also  be  defined  In  a  discrete  way.  It  was  found  that 
defining  the  discrete  operators  by 


P2  9x(rj»z) 


[^9jt(r j »z)  -  9ji(rj+Ar,z)  -  g*(rj-Ar,z)] 

(ZFp 


(6-79) 


and 


V  9i(rj,z) 


*2  -  V4 

“7T 


9jt(rj>z) 


(G-80) 


leads  to  errors  In  the  numerical  solution.  This  is  most  easily  seen  by 
considering  the  example  in  which  Jl-0  and  ao*l  so  that  9o  *  /r.  For  the 
continuous  operators  of  (G-77)  and  (G-78),  (P2+V)/r  *  0.  However,  for  the 
discrete  operators  of  (G-79)  and  (G-80) 

(Pa*V)  g,(rj)  .(^i  -  /J=T  -  /3FT  -  $1  ]  *  0  (G-81) 


This  inconsistency  is  most  pronounced  near  the  origin  where  it  may  cause 


serious  distortions. 

To  correct  for  this  problem  we  modify  the  operator  P2+V  so  that  the 
equation 


[P 2+V]  /r  JA(k nlr)  =■  (knA)2  /r  JA(knJlr)  (G-82) 

is  satisfied  exactly  when  restricted  to  the  discrete  r.-  and  the  slowest 
eigenmode  VT  JA  (k 0jj.  r).  If  we  retain  the  definition  for  P2  in  (G-79), 
this  condition  leads  to  the  following  set  of  equations  for  V(rj) 

V(fj)  /rj  Jx(k01rj)  =  (k0jt)2  /rj  J*(k0Jlrj) 

Jjtfroifl)  -  /F77T  l  (G-83) 

Ar2 


The  discrete  incremental  wave  equation  is  now  propagated  by  applying  the 
Fourier  transform  and  its  Inverse  to  the  system  of  N  equations 


gx(rjtz+Az)  -  e 


if- 


V(r 


J  ^  f~  1 


p 2 

[e  Zlc  gA(rj,z)]  -  AzvF  sA(rj>Z) 


-1  Az  2f  1-cos  (2ni/N)  1 
[«T  F  9jt(rj,2)] 


Az/r  sA(rj,z 


(G— 84) 


4.  Resonant  Particle  Dynamics  in  Two -Pi mens ions 

A  2-D  computer  code  has  been  written  to  propagate  the  laser  beam 
using  fast  Fourier  transforms.  A  number  of  simplifying  assumptions  have 
been  made  In  the  code  to  permit  us  to  quickly  obtain  useful  Information 
about  the  differences  between  the  one  and  two  dimensional  analyses. 

The  code  assumes  that  the  pump  field  Is  produced  by  a  wiggler  magnet  and 
that  there  Is  perfect  Injection  of  electrons  so  that  the  electron  trans¬ 
verse  canonical  momentum  equals  zero.  It  Is  assumed  that  the  transverse 
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component  of  the  electron  velocity  is  determined  entirely  by  the  magnetic 
field.  To  study  the  effect  of  the  applied  axial  electric  field  on  the 
resonant  phase  of  the  trapped  particles  we  also  restrict  the  electron 
dynamics  to  that  of  the  resonant  particle.  In  the  resonant  particle 
description,  only  the  trapped  electrons  contribute  to  laser  gain. 

The  effective  transverse  current  density  is  therefore 

J±  =  Ft  -pj  $j_  ( G— 85 ) 


where  Ft  is  the  fraction  of  the  electron  distribution  that  is  trapped. 
From  (G-54) 


e  1  = 


Using  (G-56)  and  (G-58) 


-► 


eA  eAp 

ymc  *  "  ymc  * 


iz  _e_  ^ 

PJ  Tmc  HP 


(G-86) 


(G-87) 


where  TR  is  the  phase  of  the  resonant  particle.  Since  the  pump  field  is 
produced  by  a  wiggler  magnet  we  can  chose  an  initial  phase  $p  *  0.  Then 
Ap  *  A £  *  Bm/km,  where  Bm  is  the  amplitude  of  the  transverse  wiggler  field 
and  km  is  the  wave  number.  The  wave  equation  (G-60)  then  becomes 


ii 


VX2 


+  i 


kl  sfl 


A,  e 


1* 


2  ymc 


e  J, 


(G-88) 


From  Maxwell  's  equations 


(G-89) 


Differentiating  (G-4)  and  making  the  slowly  varying  amplitude  and  phase 
approximations  we  obtain 


aT-  "  1  “l  *l  <6-*» 
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which  when  substituted  into  (G-89)  gives 

K  ■  fL  • 

Using  (G-91)  in  (G-88)  we  obtain 

WL7l2  -  -5z ’]  ELe  *L 

which  is  the  equation  that  is  solved  in  the  2-0  computer  code. 

In  the  code  the  laser  electric  field  is  specified  by  two  column 
vectors  of  N  elements  each,  E^rj)  an<*  \(rj)»  ^  = 

The  source  term  that  produces  laser  gain  is  given  by  the  right  hand  side 
of  (G-92)  which  is  also  a  column  vector  with  N  elements. 

To  optimize  the  gain  in  an  FEL  amplifier  an  axial  electric  field,  or 
alternatively  a  wiggler  taper,  could  be  used.  The  axial  field  (or  taper) 
performs  two  functions.  First,  It  compensates  for  the  phase  shift  of  the 
ponderomotive  wave,  produced  by  the  interaction  of  the  electromagnetic 
field  with  the  electron  beam,  and  second,  it  determines  the  value  of  the 
resonant  phase.  In  a  1-D  model  of  the  FEL  the  amplitude  and  the  phase 
of  the  laser  field  at  a  given  axial  position  are  uniquely  determined  and 
a  value  of  the  axial  electric  field,  Ez,  can  be  specified  which  both 
compensates  for  phase  shift  and  produces  a  desired  value  of  the  resonant 
phase.  In  a  two-dimensional  model  the  electromagnetic  field  can  have 
different  values  of  amplitude  and  phase  at  different  radial  positions 
and  the  amount  of  phase  shift  can  also  vary  as  a  function  of  radial 
position. 

Since  amplitude  and  phase  are  no  longer  uniquely  specified,  it  is 
only  possible  to  chose  one  particular  radius  at  which  both  the  phase 
shift  can  be  compensated  and  the  resonant  phase  specified.  At  all  other 
radial  positions  the  resonant  phase  must  be  calculated  and  cannot  be 
chosen  a  priori.  To  optimize  laser  gain  a  particular  radius,  r0pt«  is 
chosen  at  which  the  gain  Is  to  be  maximized,  and  a  condition  Is  specified 
for  optimizing  the  gain.  One  example  of  such  a  condition  is  that  the 
phase  space  area  of  the  bucket  at  rQpt  remain  constant  along  the  ampli¬ 
fier. 


For  simplicity  we  first  consider  the  case  in  which  phase  shift  can 
be  neglected.  The  axial  electric  field  in  the  amplifier  that  is  needed 
to  produce  a  particular  value  of  at  ropt  is  given  by 


e  ^L ( ropt ,z )  ^m(r opt ,z ^ 

:  Y^ropt»2^  ^z^ropt'2' 


sin  Vropfz>  (G‘93) 


It  is  assumed  in  (G-93)  that  £z (z)  is  independent  of  radial  position.  For 
the  segmented  waveguide  electrode  structure  we  are  considering  this  is  an 
excellent  assumption.  With  the  value  of  E(z)  given  in  (G-93)  the  value  of 
the  resonant  phase  at  all  other  radii,  rj,  is  given  by 

Vrj>  ■  sin-1  C-5-m.  T(rj)  ez(rj)  ,  (S-94) 


where  the  z  dependence  of  the  quantities  in  (G-94)  has  not  been  shown 
expl icitly. 

All  of  the  quantities  in  (G-94)  are  either  constants  or  have  a  slow 
radial  dependence,  except  for  EL  and  Therefore, 

sin  Vrj)  a^Jrp  (G“95) 

Now  in  the  resonant  particle  model,  laser  gain  Is  proportional  to 
sin  times  the  fraction  of  the  electrons  trapped  In  the  bucket.  If  we 
assume  the  phase  space  buckets  are  filled  with  particles,  then  the 
fraction  trapped  will  be  proportional  to  bucket  area.  Bucket  area  is 
proportional  to  the  quantity  E^^n(f^)  where  n(TR)  is  the  ratio  of  the 
area  of  a  bucket  of  resonant  phase  to  the  area  of  a  bucket  of  resonant 
phase  yR*0*  Making  these  assumptions 

dE.  l/_  ,  n(? o) 

.  E^n(IR)  sin  T„  -  (G-9S) 

If  the  resonant  phase  YR(rj)  *  90° ,  the  bucket  area  will  equal  zero  and 
there  will  be  no  gain.  From  (G-95)  this  means  that  laser  gain  will 
decrease  to  zero  for 

EL(rj)  4  E^(rpp^)  sin  ^(^gpt)  (G-97) 


This  Is  shown  In  Figure  17  In  which  a  resonant  phase  angle  of  24°  was 
chosen  at  ropt* 


If  the  phase  shift  is  taken  into  consideration  then  both  the  ampli¬ 
tude  and  phase  of  the  laser  beam  will  vary  as  a  function  of  radius  and  the 
formula  that  determines  the  axial  field  needed  to  produce  a  given  resonant 
phase  at  ropt  is  more  complicated  than  (G-93).  If  a  phase  shift 
of  A$  (ropt)  occurs  at  ropt  over  a  distance  A z,  the  effect  of  this  shift 
can  be  negated  by  accelerating  or  decelerating  the  electrons  over  the 
interval  A z  so  that  the  electron  distribution  remains  in  the  same  position 
relative  to  the  bucket.  The  desired  value  of  the  resonant 
phase  Y^(fopt)  is  then  produced  by  adding  to  the  field  required  to 
compensate  for  phase  shift  an  additional  amount  equal  to  the  field  needed 
to  produce  the  desired  resonant  phase  in  the  absence  of  phase  shift.  The 
resultant  value  of  Ez  is  then  given  by 


ymc 


IlJe 
tCX 


sin 


m 


+ 


Ymc2  r _ -2r2  M  + 

eAz  L(l+a2)  (kL+km)  Az 


(G-98) 


where  all  terms  are  evaluated  at  rQpt.  Y0  is  the  intial  value  of  the 
resonant  energy  at  the  input  end  of  the  amplifier  and 
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*  Bm 

ici^ 
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at  radii  other  than  ropt  the  phase  shift  will  in  general  be  different 
than  A<>  (ropt).  The  resonant  phase  at  r^  *  ropt  can  be  determined  by 
first  finding  the  electric  field  required  to  compensate  for  phase  shift, 
E'. 

The  electric  field  AEZ  that  remains  after  E'  is  subtracted  from  Ez 
can  then  be  substituted  for  Ez  in  (G-94)  to  obtain  the  resonant  phases  at 

rj  *  ropt* 

AE(r j)  -  Ez  -  E * (r j )  (G-100) 


Vj> 


(G— 101 ) 
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